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Abstract 

A mechanism for generating metric perturbations in inflationary models is 
considered. Long- wavelength inhomogeneities of light scalar fields in a decoupled 
sector may give rise to superhorizon fluctuations of couplings and masses in the 
low-energy effective action. Cosmological phase transitions may then occur that are 
not simultaneous in space, but occur with time lags in different Hubble patches that 
arise from the long- wavelength inhomogeneities. Here an interesting model in which 
cosmological perturbations may be created at the electroweak phase transition is 
considered. The results show that phase transitions may be a generic source of 
non-Gaussianity. 
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1 Introduction 



The energy density during inflation is dominated by the inflaton potential energy. At 
the end of inflation, the energy stored in the inflaton potential is converted into particles, 
which decay and reheat the Universe by thermalization to start the standard hot big bang 
phase. 

In this paper we consider phase transition from a phase A to another phase B, which 
are distinguished by the scaling of the energy density of a component pi. Namely, we 
consider a scaling of the energy density pi oc a""-* in phase A and pi oc a""^ in phase B, 
where nA ^ ub causes generation of the density perturbations when the phase transition 
is inhomogeneous in space. The mechanism is very general and can be applied to many 
other models in which the transition between phases of different scaling is inhomogeneous 
in space, even if the transition is not a "phases transition" in the strict meaning. 

Before discussing inhomogeneous phase transition, we review the mechanism of inho- 
mogeneous reheating to illustrate the basis of inhomogeneous scenarios. In Ref. P], it 
has been argued that in realistic models of inflation the coupling of the inflaton to matter 
can be determined by the vacuum expectation values of fields in the underlying theory. If 
those fields (in the string theory they would be moduli fields from the compactified space) 
are light during inflation, they will fluctuate leading to density perturbations through the 
inhomogeneities of the coupling constants. 

If the density perturbations created during inflation are negligible and the Universe 
after inflation is filled with particles ip of mass and decay rate < Hj, where Hj is 
the Hubble parameter during inflation, spatial inhomogeneities in may lead to density 
perturbations when the particles decay into radiation. In deriving the magnitude of the 
density perturbations arising from the inhomogeneity, it is useful to compare the energy 
density in a region to the virtual hidden radiation p^j^, which scales aj^ 



and calculate the density perturbations on a uniform pl^^ surface. Here, we assume that 

there is no energy transition between the radiation density pj^^ and other components of 

^The "virtual hidden radiation" is introduced just to keep track of the unperturbed spatially flat 
hypersurfaces. 
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the Universe. Assuming that the domination by ifj particles starts at adom = o-itdom) when 



Pdom = Pitdom) - P^itdom) - 

domination interval tdom <t<tdt 



the energy of the ip particles scales as matter in the 



p^c::^ poca ^ (1.2) 

with decay time t = tdec defined by 

Pdec = p{tdec)^TlMl (1.3) 

Outside the domination interval, we assume that the energy density scales as radiation. 
Fig.([T]) shows a schematic representation of the inhomogeneous boundary that creates 
density fiuctuations. Note that in this model the delay of the ■j/'-decay causes a delay 




Figure 1: Due to the inhomogeneity, the decay of the i}) particles does not occur 
simultaneously in space, which leads to a fluctuation of tdec and pdec/ p[jh- Thus, the 
evolution of the energy density is different in different patches, which results in density 
fluctuations. 

in the evolution of the energy density. The calculation of the density perturbation is 

•^Here we assume that the mass of the 'ip particles is a constant. Unhke the original argument in ref.[l], 
we consider a uniform Pdom and 6M^ = to simplify the argument. 
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straightforward. Considering the domination interval, we find 

(^dec \ Pdom Pdom / -. . \ 

where pdom and Mp are uniform in space, while is inhomogeneous. Using pj^^ in Eq.f ll.ip . 
we can find the energy density after the decay: 



1/3 

Ode^ . _ Pdom r 
P OC ^ Pvh- „2/2...2/zP'"h-> \^-^> 



where the ratio p/p^/j is a time-independent constant after the decay. The density per- 
turbation on a uniform pj]^ surface is thus given by 



5p 2 5T^ 



;i.6) 



p 3 

which reproduces the limit F^/if/ — > in Ref. [2]. 

Another way to generate cosmological perturbations from an inhomogeneous boundary 
is to consider an inhomogeneous end for the infiationary phase [3l HJ O [6] . For inflationary 
expansion, the equation for the number of e-foldings is 

TV^ln^, (1.7) 

where is the time when the long- wavelength inhomogeneity exits the horizon and is 
the time when inflation ends. We deflne (pN = <P{tN) and 0e = (t>{te) for the inflaton fleld 
0. Using pj]^ and repeating the calculation given above, in place of Eq. (11.41) and Eq. 
(II. Sp . we obtain 



a{tN)J pit 



N 



I (^[te) \ r AN r /i n\ 



and 

If we assume instant decay and instant thermalization after inflation, the energy density of 
the Universe after inflation scales as radiation and p/pl^ is a time-independent constant 
after inflation. Therefore, the density perturbation on a uniform p^^ surface, which is 
caused by the inhomogeneities in N, is given by 

^ = ASM, (1.10) 
P 



and recovers the conventional delta-N formula ( = 6Nu More specifically, we can calculate 
6N from the (pe inhomogeneity in the inflationary scenario using a very simple equation 
SNend — {dN / d(j)e)^'Pe- In niost inflationary scenarios, is given explicitly by (f)^ and Re- 
considering the two scenarios discussed above, the curvature perturbations created by 
the inhomogeneous boundaries are natural consequences of the inhomogeneities arising 
from long-wavelength fluctuations of light fields. In this paper, we consider inhomoge- 
neous phase transitions in which the critical temperature is not homogeneous in space. If 
the potential energy dominates during a short interval, the phase may be dubbed mini- 
inflation. Following the uniform p^h calculation discussed above, we can calculate the 
density perturbations created at the phase boundary. To calculate the density perturba- 
tions we assume (1) the beginning of the phase occurs simultaneously in space, but the 
end is inhomogeneous, (2) the transition occurs instantly just after the interval and (3) all 
the energy stored in the potential is translated into radiation. Complementary scenarios 
for more generic situations require numerical study and are highly model-dependent, thus 
they will be considered in future works. However, we consider a particularly attractive 
model, featuring the possibility of inhomogeneous phase transitions at the electroweak (or 
more generically, unification) scale that may lead to the creation of a significant level of 
non-Gaussianity. 

2 The model 

2.1 Simple model for second order phase transition 

To illustrate some typical features of finite temperature effects, we consider a real scalar 
field and a potential: 

V{<f>) = ^-lmJ02 + lA0^ (2.1) 

where Vq is tuned so that the cosmological constant vanishes at the true minimum. The 

phenomenon of high-temperature symmetry restoration can be understood by the finite- 

''See Appendix A for the definitfon of the curvature perturbation C, and the 6N formula that relates C 
to 6N. 
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temperature effective potential given by [7j 



rp4 

2^ 



dx In 



1 — exp 




x^ + 



y2 



(2.2) 



where V{(j)c) is the one-loop potential for zero-temperature with the classical field 



V(0c 



(2.3) 



where is a renormalization mass scale. At high temperatures, Vr can be expanded as 

1 



(2.4) 



which suggests that the temperature-corrected effective mass at 0c = changes sign at a 
critical temperature 



2m^ 
Ai/2- 



(2.5) 



In a more general situation, one may introduce couplings to the fields in the background 
thermal bath. If the couplings of to the fields in the background thermal bath are more 
significant than the self-coupling, a typical form of the potential with a thermal correction 
term is given by 

V = Vo+(^g^T^-\ml^<f + ..., (2.6) 

where g denotes the effective coupling of to the fields in the thermal bath. In this case, 
the critical temperature is given by 



(2.7) 



2^ 

In this section, we consider the latter case where is given by ~ 

The phase transition is second order in the model discussed above. We consider two 
distinct cases: 

1. The energy density of the Universe is dominated by the potential energy Vq during 
the interval Tdom > T > T^. The Universe is then dominated by radiation, due to 
instant decay. We assume that all the energy stored in the potential is converted 
into radiation just after the phase transition (i.e., we assume = T^ec- See also 
FiglSD. 
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2. The energy density of the Universe is still dominated by radiation at T = Tc. 
After the phase transition at T = Tc, all the potential energy is converted into 
non-relativistic particles ip that scale as matter. The interval of the radiation dom- 
ination may end at T = Tdom when p^/prad — 1, or more generically the ip particles 
may decay into radiation at T = T^ec before the domination. In this scenario, the 
inhomogeneous phase transition causes the inhomogeneities of the matter density. 
See also Ref.[T3] in which the inhomogeneities of the curvatons are generated by 
inhomogeneous preheating. 
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Figure 2: Initially, the Universe is dominated by radiation. The potential energy then 
starts to dominate at T = Taom- The domination by the potential ends at T = T^, where 
a phase transition occurs. Radiation domination starts after the phase transition. 

In the former case, calculating the density perturbation is straightforward. We assume 
that the interval of domination by the potential energy starts 8X T — TdQ^ = T(t(j^QYa). 
Considering plf^ oc as before, after the phase transition at T = Tc = T(tc) we find 
that 

P0W«(^)%;.- (2.8) 
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Just after the phase transition, the potential energy is converted into radiation. The 
energy density perturbation on a uniform p^;^ surface is thus given by 

^ = -4^ = -4^+ A (2.9) 
P Tc g 

where the curvature perturbation is given by 

lSp^_S^^_S^^S,^ (2.10) 

This result can be obtained alternatively from the delta-N formula ( = 6N. For infla- 
tionary expansion during the Vo-dominated interval[8], the number of e-foldings is given 
by 

N = \n('^), (2.11) 



T 

which leads to ( = 6N = —6Tc/Tc. In order to calculate the pure contribution from the 
inhomogeneous phase transition, we assume that all the energy stored in the potential is 
converted into radiation just after inflation. To understand the light-field potential, we 
consider a specific choice for the cr-dependent mass: 



m^(a) = (^1 + a— J , (2.12) 

where a is the light field and A is the cut-off scale of the effective action. Note that 
a conventional interaction ~ arn^a"^ cfP' / in the effective low-energy action may induce 
the (j-dependent mass. In this case, the thermal correction to the mass of the light 
field is m^(T) ^ a(mQ/A^)T^, which is supposed to be smaller than the Hubble param- 
eter ~ Maxlp™*^, Vo}/3Mp, as in the inhomogeneous reheating scenario discussed 
in Ref. pLj- If there is no significant potential other than the finite-temperature effective 
potential V{(f))Ti we find an effectively flat a potential during the symmetry restoration 
phase. Since the interaction depends on the values of the fields a and 0, the background 
field trajectories after the phase transition may be sensitive to the initial conditions and 
the non-perturbative effects of the decay process, which means that the general evalua- 
tion of the cosmological parameters after the phase transition typically requires numerical 
calculations [9]. However, the numerical study related to such a non-perturbative process 
after the phase transition is highly model-dependent and out of the scope of this paper. 
We thus assume that all the energy stored in the potential is instantly converted into 



radiation just after the phase transition, in order to single out the contribution from the 
inhomogeneous phase transition. In addition to the complexities of the decay process, 
the domain walls related to discrete symmetry breaking may cause a problem. However, 
cosmological domain walls can be made unstable and safe if a bias between the two vacua 
is induced by an effective interaction term that breaks the Z2 symmetry. Note that for 
supergravity, domain walls caused by R-symmetry are safe since the supergravity interac- 
tion creates the required bias [lOj. Therefore, for simplicity and to allow calculation of the 
model-independent contribution from the inhomogeneous phase transition, we ignore the 
domain wall problem in this paper, expecting that the walls decay instantly into radiation 
due to the bias between the two vacua. 

The latter scenario is less trivial. Let us consider the case in which the energy density of 
the Universe is dominated by radiation at T = and the non-relativistic i\} particles decay 
into radiation at T = T^ec < Tc-, as is shown in Fig. [31 We assume that all the potential 
energy is translated into il) particles at T = Tc- We introduce the ratio r = p^/p and 
consider the case in which i\} does not dominate the universe (i.e., r{tiiec) < !)• Assuming 
that the symmetry restoration phase starts at some uniform temperature T = Tr, and 
introducing p^^ as before, at t = tdec we find 

PAt)^^^]' (^]p:h, (2.13) 



Tc J \ Tficc 

where pj]^ scales like radiation. The decay temperature T^ec = Tit^ec) is determined by 

p{Uec)^TlMl (2.14) 
where we assume 5T^ = 0. Therefore, the density perturbation is given by 

6p STr 

- = -3r-^. (2.15) 
P Tc 

It is possible to consider a case in which the potential energy decays into radiation 

immediately after the phase transition (i.e., for r{tc) < 1 and Tc = Tdec)- In this case, we 
find 

p^t) oc (^^^ pl,, (2.16) 
which leads to the density perturbation given by 

^ = -4rf^. (2.17) 

P Tc 
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Figure 3: Initially, the Universe is dominated by radiation. The potential energy is 
converted into non-relativistic matter at T = Tc. Then the matter decays into radiation 
ai T — T^ec- In the picture we show a case in which the non-relativistic matter dominates 
before the decay, but it is possible to consider a case in which the decay into radiation 
occurs before domination, as is discussed in the text. 

2.2 Non-t hernial trapping 

In the simple second-order example, we consider effective couplings that depend on light 
fields. Long-wavelength inhomogeneities of the light fields may lead to an inhomoge- 
neous critical temperature 5Tc ^ 0. Here, we consider another example, in which long- 
wavelength inhomogeneities of the number density of some particles cause an inhomoge- 
neous end of the symmetry restoration phase. 

During preheating, some of the kinetic energy of the infiaton is converted into excita- 
tions of the preheat field If x couples to a field (pa with a potential 

V{<Pa, x)-Vo- \m'<Pl + - ^0^x^ (2.18) 
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where the inflaton terms are omitted, the effective potential caused by the high density 
of the preheat field is given byjSj Hi] 

1 

V^^f{ct>a) ^V,- + + g\^Mx- (2-19) 

For > 0, the effective potential gives 



f'"(«-V„-i™=(^J-^)VC=f. (2.20) 



When is very large, the field (pa is trapped by a strong attraction from the origin. 
During the interval of the trapping, the potential barrier decreases, since scales as 
n-^ oc a~^, and ultimately tunneling occurs below the critical number density [SlITl] given 

by 

TTl 

ny<nc = — . (2.21) 

If the energy density is dominated by the potential energy Vq, the trapping leads to an 
inflationary expansion. The number of e-foldings elapsed during this interval is 

N = Un(^], (2.22) 



3 \n^{te 

where ti and te are the time when the domination by the potential energy starts and when 
the inflationary expansion ends. In this model, inhomogeneities in the initial number 
density n^iU) can be created by inhomogeneous preheating [121 [T3]lf| The inhomogeneities 
in the preheating arise from the long- wavelength fluctuations of the multi-field trajectory 
for the symmetry-breaking potential. This is the origin of 6N discussed in Ref.|i5j. In 
addition to the inhomogeneities in 6n^{ti), we may consider inhomogeneities in n^, which 
are non-zero if m and g are modulated at the end of the trapping phase. Using the delta-N 
formalism, for Sn^{te) — Suc we find 

5iV = -i^^-^ + l^. (2.23) 
3 Tie m 3 g 



^Inhomogeneous preheating accompanied by instant decay may directly lead to the creation of cur- 
vature perturbation|12] . In this section, we consider a preheating field that does not lead to instant 



decay [IS US]. 



11 



2.3 Electroweak phase transition 

The main obstacle in building a model of an inhomogeneous electroweak phase transition 
is that the long-wavelength inhomogeneities of the light field must survive until the elec- 
troweak phase transition, when the Hubble parameter is much lower than the gravitino 
mass. In supergravity models inspired by string theory, there are many light fields (mod- 
uli) in the effective action, but typically the mass of the moduli fields is expected to be 
of the same order as the gravitino mass, where the gravitino mass is generically given by 
m^l2 ~ Ag^gy/Mp, where Asusy > TeV is the supersymmetry breaking scaleP4j. The 
mass of the moduli ~ ^3/2 is clearly larger than Hew = T'^wl^p^ where Tew is the 
critical temperature for the electroweak phase transition. Therefore, if an inhomogeneous 
phase transition occurs at the electroweak phase transition, the inhomogeneities of the 
effective action must be inherited from the fluctuations of the light flelds whose poten- 
tial is protected by symmetry before the electroweak phase transition, while the moduli 
potential must be lifted after the phase transition. The above condition for the inhomo- 
geneous electroweak phase transition might seem very severe, but in string theory there 
is at least one speciflc example that may induce an inhomogeneous phase transition at 
the electroweak scale. We consider intersecting D-brane models, which are an interest- 
ing possibility for string model building, allowing us to devise models that are sensibly 
close to the Minimal Supersymmetric Standard Model (MSSM) in terms of particles and 
gauge groups [TSj. A remarkable feature of this scenario is that the flavor structure of the 
Yukawa couplings may arise from the matter flelds located at different intersections, with 
the resulting Yukawa couplings expressed by the classical instanton action of the minimal 
world- sheet 



where A is the minimal world-sheet area of the intersection. If the model is constructed 
from D6-branes in Type IIA string theory wrapping orientifolds of i?4 x x x T^|15j. 
there will be shift symmetries that correspond to the brane motion in the internal space. 
If the shift symmetries are not broken in the effective action, the minimal world-sheet 
area remains as an arbitrary parameter. Considering moduli fields cTj, [i = 1, 2, 3) for the 




(2.24) 
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three branes constituting a triangle in the internal space, we find; 



(2.25) 



da id (J j 



It would be better for our purpose to consider a simple form of A{cr) and consider the 
inhomogeneity 6a to obtain 



Let us consider a possible mechanism for generating an effective potential related to A. 
Assuming that the Yukawa couplings are generated by the mechanism and considering 
the standard one-loop correction to the Higgs field potential from the top fermion loop, 
we obtain for 



where m^^ denotes the scalar top mass[T7j. Here we consider Yf oc exp(— A/(27ra')). In 
the MSSM, the one-loop correction from the top Yukawa coupling destabilizes the Higgs 
potential and causes electroweak symmetry breaking. From the one-loop correction term 
in Eq. ( \2.27\\ . we find that the world-sheet area can be stabilized after the electroweak 
symmetry breaking p!6]. In this case, the free motion of the D6-branes in the internal space 
is protected by the shift symmetry before the electroweak symmetry breaking. However, 
after the electroweak symmetry breaking, which is induced by the loop correction in the 
MSSM electroweak symmetry-breaking scenario, the shift symmetries are partly broken 
and the minimal world-sheet area is stabilized in the low-energy effective action. Although 
the scenario depends greatly on the specific details of the intersecting brane models, a 
generic implication of the scenario is that the inhomogeneous phase transition occurs 
whenever the shift symmetries are not explicitly broken before the phase transition. A 
similar mechanism may work at the GUT phase transition, and the inhomogeneous phase 
transition may lead to a cosmological signature of the intersecting brane models. 

3 Conclusions and discussions 

We have studied a mechanism for generating primordial density perturbations in infia- 
tionary models. We considered long-wavelength inhomogeneities of light scalar fields that 




(2.26) 




(2.27) 
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cause superhorizon fluctuations of couplings and masses in the effective low-energy ac- 
tion. Since the effective couplings and masses are not homogeneous in space, cosmological 
phase transitions may occur that are not simultaneous in space. It is possible to create 
the primordial curvature perturbation from the mechanism, but more generally the sce- 
nario of an inhomogeneous phase transition allows for non-Gaussianity to occur in the 
spectrum after inflation [HI [19]. It is useful to specify the level of non-Gaussianity by the 
non-linear parameter /nl, which is usually defined by the Bardeen potential $, 

$ = ^Gaussian + InL^ Gaussian' (3-1) 

Using the Bardeen potential, the curvature perturbation ( is given by 

$ = ^C- (3.2) 

When we consider "additional" non-Gaussianity created at the inhomogeneous phase 
transition, the first-order perturbation is generated dominantly by the usual inflaton per- 
turbation. Therefore, the "additional" second-order perturbation is not correlated to the 
first-order perturbation. In this case, the non-linear parameter is estimated as [20] 



IfNL ^ ^ K'^- + NLV^ \ogihL)] , (3.3) 

where ( can be expanded by the 6N formalism as 

C ~ N^6<f) + NJa + ^N^^6<j)^ + ^N^Ja^ + (3.4) 

and we assume that the perturbation can be separated as 

^ = + (3.5) 

Here fcf, = min{A;j} {i = 1, 2, 3) is the minimum wavevector of the bispectrum and L is 
the size of a box in which the perturbation is defined. A useful simplification is|21j 

The scenario of adding non-Gaussianity from the inhomogeneous phase transition is in- 
teresting, since for the effective low-energy action, higher-dimensional couplings may nat- 
urally appear with light fields in a decoupled sector. 
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Consider a simple example discussed in Sec. 2.1 with (5m^ 7^ and 6g = SX = 0. 
Considering the initial value for the light field a in Eg. (12. 121) . a modest assumption would 
be cr ~ 0. From Eq. fl2.10p . the curvature perturbation created from the inhomogeneous 
phase transition is purely second order and given by 

A.) ^ ^ = (3 7) 

^ ~ " 2A2 2A2(27r)2' ^^''^ 

where Hj is the Hubble parameter when the long-wavelength inhomogeneity of the light 
field a exits the horizon during infiation. Thus we find from the 6N formula; 

iVa. = (3.8) 

Even for the initial condition a ^ 0, the non-linear parameter for the inhomogeneous 
phase transition is significant. Considering the usual normalization for the first order 
perturbation, we find[S] 

\N^6(f)\ ^ 5 X 10"^ (3.9) 
The non-linear parameter is thus given by 

fNL^ (^lO^xall^ . (3.10) 

Considering the modest bound for the non-linear parameter I/nlI < 100, the above result 
puts a significant upper bound on the infiationary scale or on the effective couplings that 
contain decoupled light fields. 

For the electroweak phase transition, we find for the simple case {A = A{(t)); 

^ AruH 2Tia' 2Txa' Ka 2TTa' \Aa J ' 

With regard to the non-Gaussianity, we find from the above equation that a in the 
standard calculation is simply replaced by — a2/(2vra') for the electroweak phase transition 
with the effective scale = A. 

Inhomogeneous phase transition can be applied to warm inflationary scenario. In- 
homogeneous phase transition in warm inflation will be discussed in our forthcoming 
paper [22]. 
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A SN formalism for the curvature perturbation 

Here we consider two different definitions for the curvature perturbations. The comoving 
curvature perturbation (TZ) can be related to the curvature perturbation on uniform- 
density hypersurfaces {() by studying the evolution at large scales. The gauge-invariant 
combinations for the curvature perturbations can be constructed as follows: 



(j) with a standard kinetic term. Linear scalar perturbations of a Friedman-Robertson- 
Walker(FRW) background were considered: 

ds^ = -(1 + 2A)dt^ + 2a\t)ViBdx'dt + a\t)[{l - 2ij)-fij + 2V iV jE]dx'dxK (A.2) 

Here p and p denote the energy density and the pressure. Spatially flat hypersurfaces and 
uniform density hypersurfaces are defined by = and 5p = 0, respectively!^ 

Besides the curvature perturbations defined above, it is useful to define the perturbed 
expansion rate with respect to the coordinate time. The perturbed expansion rate is 
expressed as 



Choosing the gauge whose slicing is fiat at tini and uniform density at t, the 6N formula 
is given by 




56 = -3^ + V V, 



(A.3) 



where the scalar describing the shear is 



a = E-B. 



(A.4) 




(A.5) 



^In this Appendix, "-0" is used for a metric perturbation. 



16 



The 5N formula is sometimes expressed by 



C = 6N = -H^-P 
P 



(A.6) 



where 6N is the perturbed expansion to uniform- density hypersurfaces with respect to 
spatially flat hypersurfaces, and 6p must be evaluated on spatially flat hypersurfaces. 
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